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We study the t-J-V model beyond mean field level at finite doping on the triangular lattice. The 
Coulomb repulsion V between nearest neighbors brings the system to a charge ordered state for V 
larger than a critical value V c . One-particle spectral properties as self-energy, spectral functions and 
the quasiparticle weight are studied near and far from the charge ordered phase. When the system 
approaches the charge ordered state, charge fluctuations become soft and they strongly influence 
the system leading to incoherent one-particle excitations. Possible implications for cobaltates are 
given. 
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I. INTRODUCTION 

Strongly correlated electronic systems continue to be 
at the center of interest in condensed matter, the t-J 
model being a paradigm for theoretical studies. Since the 
early days of high-T c superconductivity the t-J model on 
the triangular lattice attracted a great deal of attention, 
since it was believed to be a candidate for realizing the 
resonating valence bond scenario^. 

Recently, the interest about electronic correlations on a 
two dimensional triangular lattice received a new motiva- 
tion. Since superconductivity was discovered in hydrated 
cobaltates^ (Na x Co02-yH20) an enormous amount of 
attention has been focused on this system. Cobaltates 
are 3d-electron systems having a quasi-two-dimensional 
structure of C0O2 layers, where the Co-ions are located 
on a triangular lattice. The interplay between electronic 
correlations and the proximity to a charge density wave 
were proposed recently as important ingredients for de- 
scribing superconductivity in such materials 3 -^^. This 
proposal was motivated by various experimental reports 
which suggest the proximity of the system to charge 
orderingi^. A minimal model for strongly correlated 
electronic systems close to charge ordering is the t-J-V 
model on the triangular lattice, where V is the Coulomb 
interaction between nearest neighbor sites. The parame- 
ter V is responsible for triggering the charge order insta- 
bility. 

In this paper we concentrate on the consequences of 
the proximity to charge ordering on one-particle spec- 
tral properties. For this study we use a recently devel- 
oped large-A approac h 10 : 11 that leads to a systematic 
treatment of fluctuations beyond the mean-field level, al- 
lowing thus the calculation of the self-energy. Related 
expansions were used in the past in the frame of a slave 
boson formulation for the t-J models. However, due to 
the gauge-fields inherent to slave bosons, such an expan- 
sion cannot be performed in a systematic way 13 . Our 
formalism is based on A-operators that do not introduce 



gauge-fields, such that the present treatment is free of 
the above mentioned difficulties. In general, large N 
expansions give preference to one kind of fluctuations 
over others. In our case, the 1/N expansion (with the 
same N = 00 limit as in slave bosons) gives preference 
to charge fluctuations over magnetic ones. Therefore, it 
fails to describe the limit of zero doping of an antifer- 
romagnet on a square lattice, where antiferromagentic 
long-range order should be present, a difficulty that is 
shared also by the U(l) slave boson formulation 13 . How- 
ever, in the present case we focus on rather high doping 
(~ 30%) on a frustrated lattice, such that in fact, charge 
fluctuations arc the dominant factor. For the square lat- 
tice, as discussed in Ref.|lll| our method is better for 
large doping, corresponding to the overdoped region of 
high-T c cuprates. 



Our results show how, starting from an already re- 
duced value by correlation effects, the quasiparticle 
weight vanishes on approaching the charge ordering in- 
stability, leading to a redistribution of the one-particle 
spectral weight. Furthermore, the soft charge modes re- 
sponsible for the phenomena above are identified. 



The paper is organized as follows. In Sec. II, a sum- 
mary of the formalism and details of the self-energy calcu- 
lation are given. In Sec. Ill, the charge instability is stud- 
ied. In Sec. IV, results on one-particle spectral properties 
(self-energy, quasiparticle weight, spectral functions) are 
presented. In Sec. V, possible implications for cobaltates 
are discussed. Section VI presents conclusion and discus- 
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II. THE MODEL AND THE THEORETICAL 
FRAMEWORK 



<ij>,a 

+ jY d {S i S j --n i n j ) + vY J 



(1) 
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where t, J and V are the hopping, the exchange inter- 
action and the Coulomb repulsion, respectively, between 
nearest-neighbor sites denoted by (ij). c\ a and Ci a are 
the fermionic creation and destruction operators, respec- 
tively, under the constraint that double occupancy is ex- 
cluded, and rii is the corresponding density operator at 
site i. 

As described in Ref. [6[ the Hamiltonian {1} can be 
written in terms of Hubbard operators 14 as: 

t 



H = ~ £ (Xfxr + h.c.) 



<ij>,p 



<ij>;pp' 



v ^ xrx^-^x 

<ij>;pp' i,] 



VP 
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where, in addition, in order to perform a large- TV expan- 
sion, the spin index a was extended to a new index p 
running from 1 to N. In order to obtain a finite theory 
in the iV-infinite limit, we rescaled t, J and V as t/N, 
J/N and V/N, respectively. We included in eq. @ the 
chemical potential /i. The operators X pp are boson-like 
while X 0p and X p0 are fermion-like 1 ^. 

In Ref. Q the large- N formalism was described in de- 
tail and here a short summary is presented which will 
be useful for the calculation of the self-energy and the 
spectral function. The Feynman rules of the method are 
summarized in Fig. [TJ 

To leading order in 1/N, we associate with the N- 
component fermion field f p a propagator connecting two 
generic components p and p' (solid line in Fig[T]) 
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k and oj n are the momentum and the fermionic Mat- 
subara frequency of the fermionic field, respectively. The 
fermion variables fi p are proportional to the X-operators 



The t-J-V model on the triangular lattice is given by (/, 



X® p ) and are not associated with the 



%p ~ VNn>' 
spinous from the slave boson approach. 
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FIG. 1: Summary of the Feynman rules, a) Solid lines repre- 
sent the propagator (eq. @)- Dashed lines represent the 
8x8 boson propagator (eq. (|6])) for the 8-component 

field 5X a . Ag p and A^jJ represent the interaction between 
two fermions f p and one or two bosons SX a respectively, b) 
Irreducible boson self-energy II a 6. Double dashed lines cor- 
respond to a dressed boson propagators, c) Contributions to 
the electron self-energy E(k, uj) in 0(1/N). 



The mean-field values ro and A must be determined 
by minimizing the leading order theory. From the com- 
pleteness condition 6 , X 00 + ^ p X pp — A, r is equal to 
x/2 where x is the doping away from half- filling. On the 
other hand, the expression for A is 



J 1 



2N S 3 <H 



^ cos(k n )n F (E k ) 



(5) 



where np is the Fermi function and N s is the number of 
lattice sites. For a given doping x; [i and A must be deter- 
mined self-consistently from (1—x) = -S~Y2k n F{Ek) and 
eq. ([5]). k v is the proyection of k over the different bond 

directions 771 = (1,0), 772 = (|, ^) and 773 = (-±, ^) of 
the triangular lattice. 

We associate with the eight component boson field 
6X a = {5R , 5X, r r '\ , r'? 3 , A V1 , A m , A m ) , the 
inverse of the propagator, connecting two generic com- 
ponents a and b (dashed line in Fig. [I]), 
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where 7 q = (2V — J)(x/2) 2 ^^cos/c^ and the indices 
a, b run from 1 to 8. q and v n are the momentum and 
the Bose Matsubara frequency of the boson field, respec- 
tively. 

The first component 8R of the SX a field is connected 
with charge fluctuations via A 2 00 = Nr o(l 4- SRi), where 
X®° is the Hubbard operator associated with the number 
of vacancies at site i. 8\'\s the fluctuation of the Lagrange 
multiplier Xi associated with the completness condition. 
r/ and correspond, respectively, to the amplitude and 
the phase fluctuation of the bond variable = A (1 + 

Fermions f p interact with the boson SX a via three and 
four leg vertices, namely A pp and K pp b respectively. The 
explicit expressions for these vertices can be found in Rcf. 
0- 

Each vertex conserves momentum and energy and they 
are of 0(1). In each diagram there is a minus sign for 
each fermion loop and a topological factor. 

The bare boson propagator -D(o) a & (the inverse of eq. 
©) is 0{l/N). From the Dyson equation, D~£ = 
D^ ab — n a (,, the dressed components D a b (double dashed 
line in Fig.nj)) of the boson propagator can be found after 



the evaluation of the 8x8 boson self-energy matrix II a fc. 
Using the Feynman rules, Ii a b can be evaluated through 
the diagrams of FigQJ) as shown in Ref. [||. 

This formalism is used here for calculating self-energies 
and one-particle spectral functionsii. The Green's func- 
tion J3J corresponds to the iV-infinite propagator which 
includes no dynamical corrections; these appear at higher 
order in the 1/N expansion. Using the Feynman rules, 
the total self-energy in 0(1/N) is obtained adding the 
contribution of the two diagrams shown in Fig. [TJ:. 

After performing the Matsubara sum and the analyt- 
ical continuation ioj n = lu + i£, the imaginary part of S 
is 

im£(k,w) = -^-^2h a (k,q,uj - E k _ q ) 
s q 

5 afc (q, uj - E k _ q ) h b (k, q,0J- E k _ q ) 
x[n F (-E k - q ) + n B {uj - E k - q )] (7) 

where n B is the Bose factor, B ab is 

B ab (q,v) = -2 1im Im [D ab (q, w n -> v + irj)}. (8) 

"q— >0 

and the 8-component vector— h a (k, q, v) is 



h a (k,q,u) 



2E k - q - 



v 



- q -f);-2A cos(k m - q -f); 



} ^ + 2A^cos(A ;n -|)cos(|);l;-2A co S (k m 
-2A cos(fc, 3 - ^);2A sm(k m - ^);2A sin(fc 1)2 - ^);2A sm(k m - ^)) . 



(9) 



It is interesting to show Im £ expliciting the terms 
Brr, B\r and B\\\ 

Im E(k,w) - {n 2 B RR (q 7 u; - E k ^ q ) 

s q 

+ 2 n B XR {q,iu - E k _ q ) + B xx {q,oj- E k _ q )} 
x [n F {-E k - q ) + u b (uj - E k - q )] + /(k,w) , 

(10) 

where fl — (E k - q + uj + 2fi)/2 + 2AV cos(fc,, — 



q v /2) cos(q v /2). In the function /(k, to) the other terms 
of eq. ([7]) have been included. After performing numer- 
ical calculations for the parameters used in the present 
paper, /(k, uj) results at least an order of magnitud lower 
than the other terms of eq. (jTOJ) . 



Using the Kramers-Kronig relations, i?e£(k,w) can 
be determined from /mS(k, w), eq. |7J and the spec- 
tral function A(k,u>) = — —Im G(k, oj) can be computed 



A(k, w ) = -I 



JmE(k,w) 
[w - £ k - fie £(k, w)] 2 + [Im £(k, w)] 2 



(11) 



Notice that the self-energy is calculated using the prop- 
agator G(k, oj) for the X-operators, that are proportional 
to the fermionic /-operators, which cannot be related to 
usual fcrmions. 



III. CHARGE DENSITY WAVE INSTABILITIES 

The inclusion of a nearest neighbors Coulomb repulsion 
V favors a charge density wave (CDW) state. In this 
section we will discuss the charge instability 6 - which will 
be of interest for analysing self-energy corrections. 

The charge-charge correlation function Xij can be writ- 
ten as±&±£ 



X 



^)4e<^w(°)> 



(12) 



Using the completeness condition and the relation be- 
tween X?° and 8R % (A 00 = Nr (l + SRi)), in Fourier 
space, 



X c (q,u;) = -7V(!) 2 D RR (q,Lu) 



(13) 



Then, in our formulation, the charge correlation function 
is proportional to the component (1,1) (also called Drr) 
of the dressed boson propagator D^. 

The divergence of the static charge susceptibility (fl3|) 
marks the onset for the CDW instability. For the tri- 
angular lattice the charge susceptibility diverges at Q = 
(4/3tt, 0) for V = V c . This new phase is called V3 x %/3 
CDW 3 ' 5 ' 6 . In what follows all energies (V, J, etc.) are 
given in units of t. 

Figure [5] shows the phase diagram in the V c — x plane at 
temperature T = 0. As in Ref. [g], we will use J = 0.2. 
The different phases are the homogeneous Fermi liquid 
(HFL), the V3 x V3 CDW, the bond-order phase (BOP), 
and a phase separation (PS) region. However, in this pa- 
per we focus only on the proximity to the HFL-CDW 
transition (see Ref. ® for details about the BOP and PS 
phases). For the large doping studied here results are 
very robust against different J values as long as they are 
not considered to be unphysically large. Only for unphys- 
ical large values J ~ 1, the BOP and PS regions approach 
the doping levels considered here. But for J < 0.5, a 
value well beyond the physical ones in e.g. cobaltates, the 
HFL and CDW phases are well defined beyond x ~ 0.15. 

Next, we study the charge dynamics on approaching 
the x 73 CDW from the HFL. In Figs. [| and H re- 
sults for im x c (Q,cj) are presented for two commensu- 
rate doping values x = 2/3 and x = 1/3, respectively. 
The momentum was fixed at Q = (4/3tt, 0). 




0.1 0.2 0.3 0.4 0.5 0.6 
doping x 



FIG. 2: Phase diagram in the V c — x plane for J = 0.2. The 
homogeneous Fermi liquid (HFL) is unstable against phase 
separation (PS), bond order phase (BOP) and the y/S x Vo 
CDW. BOP and PS are discussed in Ref. |6j and are presented 
here only for a complete description of the phase diagram. For 
a given doping x the static charge susceptibility diverges at 
V c where the onset to a y/3 x \/3 CDW occurs. 



For x = 2/3 (Fig. [3} and for V = (solid line), we 
clearly see a collective peak at ui ~ 3.lt at the top of 
the particle-hole continuum (inset). With increasing V 
the collective peak becomes soft and accumulates spec- 
tral weight as shown for V — 2.0 (dotted line), V = 3.0 
(dashed line) and V — 3.2 (dotted dashed line). At 
V = V c ~ 3.6, the collective peak reaches zero fre- 
quency condensing a V3 x v3 CDW phase. Therefore, 
the charge instability can clearly be seen as the soften- 
ing of the collective charge mode. For x = 2/3, charge 
dynamics is very similar to the case of the square lattice 
at one quarter filling where the softening of the charge 
mode was obtained and discussed in the context of or- 
ganic materials^. 

For x = 1/3 (Fig. 2]), the situation is somewhat dif- 
ferent. For V = (solid line) the collective peak is lo- 
cated at lo 2.2t, above the particle-hole continuum. 
With increasing V towards V c ~ 2.35, a transfer of spec- 
tral weight takes place over a large energy range since 
the collective charge mode spreads over the particle-hole 
continuum. 

The reason for the different behavior between x = 1/3 
and x = 2/3 is due to the form of the particle-hole con- 
tinuum for each doping (insets in Figs. [3] and [5]). For 
x — 2/3 there is a gap at low energy, since for x — 2/3 
2kp < Q, and particle-hole transitions are not possible 
for lo < 1. For V = the collective peak is formed at 
the top of this continuum (Fig. [3] solid line). When V in- 
creases, the collective peak softens (Fig. [3]) and for large 
V, near V c , emerges from the continuum until it reaches 
u> = exactly at V c (see for instance solid and dashed 
lines in Fig. [3]). In contrast, for x = 1/3 there is no 
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FIG. 3: (color online) Im,Xc(Q,ui) for x = 2/3, J = 0.2 and 
for different values of V approaching V c ~ 3.6. Q = (4n/3, 0), 
where the static charge susceptibility diverges. For V — 0, 
the collective charge peak is formed at the top of the particle- 
hole continuum (inset). With increasing V the collective peak 
becomes soft and emerges from the bottom of the continuum 
accumulating spectral weight. For V — V c ~ 3.6 this soft 
collective peak reaches u = freezing the x \/3 CDW 




CO/t 



FIG. 4: (color online) Similar to Fig. 0] for as = 1/3 (V c ~ 
2.35). Notice the different behaviour with respect to x = 2/3 
(see text). Inset: particle-hole continuum. 



gap in the particle- hole continuum (inset Fig. 2} . When 
the collective mode enters the particle-hole continuum, a 
spread of spectral weight over a wide range takes place, 
and for V approaching V c , it emerges from the continuum 
as a resonance. 

Before closing the section, we present in Fig. [5] the 
T-V phase diagram for x = 1/3 and x — 2/3. Tco is 
the temperature for charge ordering. The obtained reen- 
trant behavior of the charge order transition was also 
predicted by several calculation a 19 ' 20 i 21 ' 22 ' 23 . Recently, 




FIG. 5: T-V phase diagram for x = 1/3 (solid line) and 
x — 2/3 (dashed line). 



it was proposed that a reentrant behavior near charge 
order is necessary for describing anomalous optical fea- 
tures in 1/4- filling organic materials 2 ^ (see Sec. V for 
more details). 

In the next section, we will study self-energy correc- 
tions approaching the y3 x y3 CDW phase and discuss 
the role played by the soft modes mentioned above. 



IV. ONE-PARTICLE PROPERTIES BEYOND 
MEAN-FIELD LEVEL 

A. Quasiparticle weight 

The self-energy allows the calculation of the quasipar- 
ticle (QP) weight Z = 1/(1 - §§) at the Fermi level - 
In Fig. O the QP weight is plotted as a function of V 
for x = 1/3 (solid line) and for x = 2/3 (dashed line). 
The Fermi vectors are located in the T — K direction of 
the Brillouin zone for both, x = 1/3 and x = 2/3. As we 
will discuss below, our self-energy is very isotropic on the 
Fermi surface (FS) such that Fig. [S] is representative for 
all points on the FS. In both cases, when V approaches 
the corresponding V c , the QP weight decreases suggest- 
ing that it tends to zero at V c (see dotted arrows). Below, 
the analysis of spectral weight transfer at small energy 
will give additional support for this result. 

The behavior presented in Fig. [6] indicates the break- 
down of the Fermi liquid (FL) at V c . Therefore, our 
calculation suggests a metal insulator transition where a 
rather isotropic gap should exists for V > V c . Our ob- 
servation is based on the fact that for V < V c the QP 
weight Z is very isotropic on the FS such that the FL 
breaks down for all points on the FS. 

In both curves, when V increases the QP weight de- 
creases, first relatively slowly, and then drops to zero 
close to V c . Therefore, near the CDW instability the QP 
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FIG. 6: The quasiparticle weight Z for J = 0.2 evaluated on 
the Fermi surface as a function of V for x = 1/3 (solid line) 
and x = 2/3 (dashed line). With increasing V both curves 
show that, first, Z decreases slowly and, it drops to zero close 
to V c (Z — » for V — > V c ). Dotted lines are guides for the eyes 
which show the mentioned tendency (see text for discussions) . 
For V = the quasiparticle weight for x = 1/3 is smaller 
than for x = 2/3 which means that correlations are stronger 
for small than for large x. Inset: The QP weight divided by 
the QP weight at V = (Z/Zo) versus V/V c for x = 1/3 
(solid line) and a; = 2/3 (dashed line). With increasing V, 
the decreasing of Z is somewhat faster for x = 2/3 than for 
a; = 1/3. See text for discussions. 



carries very small weight making the electronic dynam- 
ics very incoherent. To clarify this point better we will 
present below results for the spectral functions. 

The QP weight for V = (far from the CDW instabil- 
ity) is already reduced from one. This reduction is not 
related to the CDW instability but it is due to electronic 
correlation effects which our method is able to capture 
in the pure t-J model. In fact, Fig. [6] shows that the QP 
weight for x = 1/3 is smaller than for x = 2/3, as one ex- 
pects, since correlation effects should be weaker for a di- 
lute system. For x — ► 1 our method predicts Z — > 1 which 
corresponds to the uncorrelated limit (empty band). On 
the other hand, in both cases, Z strongly decreases on 
approaching the charge order intstability. 



B. Self-energy 

The discussion above indicates that the HFL is 
strongly modified near the CDW instability, where the 
carriers are renormalized by interacting with soft charge 
fluctuations. In order to discuss this point in detail, we 
will present results for the self-energy E(k, ui). 

In Fig. [3 self-energy results at k = kp are presented 
for x = 1/3 and several values of V approaching V c . 
Panel (a) shows — TmE as a function of frequency. As 
discussed in Ref. [ll|, for V = 0, TmE is strongly asym- 




FIG. 7: (color online) Self-energy E(kF,tj) for x = 1/3 and 
J = 0.2 for several values of V approaching V c ~ 2.35. (a) 
-lmS(kF,w) versus u. It presents large structures at ener- 
gies of the order of t and shows a ~ u> 2 behavior around lo — 0. 
With increasing V a redistribution of the weight takes place, 
(b) Re E(kF,w) versus lo for the same values of V as in panel 
(a). Inset shows that, with increasing V, the slope of ReT, 
near u> = increases leading to a decrease of the quasiparticle 
weight. 



metric with respect to uj — and behaves as ~ to 2 at 
small co's. The self-energy presents large contributions 
at large energy of the order of t below the Fermi energy. 
With increasing V a redistribution of the weight takes 
place and 7m E develops structures at energies close to 
but above the Fermi energy. As Im S and Re E are re- 
lated to each other by a Kramers Kronig transformation, 
the changes in 7m E are reflected in 7?eE (panel (b)). 
For instance the slope of Re E at u> — increases with 
increasing V leading to a decrease of the QP weight dis- 
cussed in Fig. [SI The inset in panel (b) shows this be- 
havior in a smaller scale near oj — 0. While for V < 2.2 
a gradual change in the slope is observed, very close to 
V c , a stronger increase is observed, as shown by the thick 
solid line corresponding to V — 2.34. 
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FIG. 8: (color online) Spectral lunction A(k.F,Lo) versus ui for 
J = 0.2 and several values of V approaching V c ~ 2.35. For 
V = (solid line), pure t-J model correlation effects leads 
to a quasiparticle peak with reduced weight Z ~ 0.56 and 
incoherent structures at positive and negative energies. With 
increasing V, while the QP weight decreases (Fig. |BJ, weight 
is transfered to the range lo ~ 0.5t — 2.5t. 



C. Spectral functions and total density of states 

Using Im £ and Re E, the spectral function A(k, w) 
can be calculated as usual. Figure [5] presents results for 
A(\tp, lo) for several values of V. The peak at u = cor- 
responds to the QP while the other features, for instance 
at lo ~ — 3i and cj ~ 2.5t, are of incoherent character. 
For V = (solid line) the QP weight is Z ~ 0.56 (Fig. 
[6]) and the incoherent structure is mainly concentrated 
in the pronounced feature at lo ~ — 3.5i. Increasing V 
up to y = 2.2 (dotted dashed line) the QP weight de- 
creases [Z ~ 0.42) and, at the same time, spectral weight 
appears at lo ~ 0.5t — 2.5t (Fig. [8]). In addition, with in- 
creasing V, the incoherent structure at lo ~ — 3.5t moves 
to the right while loses some weight. 

In Fig. [5] it was assumed, supported by the Luttinger 
theorem, that kp does not change with the interaction V. 
On the other hand, if we wish to compare with ARPES 
experiments (see Sec.|V]), our input is a tight binding dis- 
persion which reproduces the measured FS which already 
contains the interactions. The comparison between our 
approach and Lanczos diagonalization presented in Ref. 
[Ill ] gives an additional support fort this assumption (see 
also Ref. [Hj]). 



D. Physical origin of the self-energy 
renormalizations: a 2 F(u>) 

In the following we discuss the interaction of the 
bosonic excitations described by Drr (see eq. (IT5|) ) with 
electrons that lead to the self-energy renormalizations 
discussed in the previous section. In many body theory 24 



FIG. 9: (color online) a 2 F(u>) for the same values of V as 
in Fig. [7J The features marked with solid arrows are related 
with collective charge fluctuations which become soft with in- 
creasing V and, at the same time, a 2 F(u>) acquires structures 
at low energies (dashed arrow) which lead to the decreasing 
of Z (see text). 

the quantity which contains this information is o?F(lS) 
where, F(u>) is the density of states of a boson which 
interacts with the fermions with strength a. At T = 02£, 

JmE(w) = / dn a 2 F(n), (14) 
Jo 

i.e. a 2 F(Lo) = dIm'E(uj)/du>. In eq. fTJJ the average 
over all momentum transfer q = ki — k2, where ki and 
k2 are two Fermi vectors, is assumed 24 . 

In Fig. [5] we have plotted a 2 F for the same values of 

V as in Fig. [7] As we mentioned above, our self-energy 
is very isotropic on the FS such that o?F evaluated at 
kp is representative for the average over the FS. 

Let us discuss first the case for V = (solid line). We 
can see two structures, at lo ~ — 2t and lo ~ 2t (solid 
arrows) which are related with collective charge excita- 
tions (see below). The features at lo ~ —St and lo ~ 2.5t 
(dotted arrows) are due to the fact that Im S (Fig. [7]) 
is concentrated in a finite range of energy (between the 
borders ~ — At and ~ it) so that, when performing the 
derivative dim Ti(u>)/du to obtain a 2 F{u>), these two fea- 
tures are created at these two borders. In Fig. [TO] we have 
plotted the average of the Jmx c (q, w) for all q such that 
q = ki — k2 where ki and k2 are two Fermi vectors. For 

V = (solid line) , the well defined peaks (solid arrows) at 
lo ~ — 2t and 2t are due to collective charge excitations 
while tails are due to the particle-hole continuum (see 
solid line in Fig. 2]). Clearly, the two structures marked 
with solid arrows in Fig. [S] are correlated with the cor- 
responding ones in Fig. [TO] showing that these structures 
in a 2 F are related to collective charge fluctuations while 
tails, are due to the particle-hole continuum. The rea- 
son why a 2 F(u}) does not trace exactly the average of 
the charge correlations of Fig. [TO] is the following. a 2 F 
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FIG. 10: (color online) Average of lm^ c (q,a)) for all q such 
that q = ki — k2 with ki and ki on the FS. With increasing 
V, the collective charge peaks (solid arrows) become soft and 
at the same time a soft dynamics (dashed arrows) appears 
at low energies. The close correlation between these features 
and those discussed in Fig. [9] shows that the charge dynamics 
is the main cause for the reduction of Z near charge order. 



FIG. 11: (color online) Spectral functions at x = 1/3 in the 
main directions of the Brillouin zone for (a) V — and (b) 
V = 2.2. The peak near u — is the coherent quasiparticle 
peak and the other structures are of incoherent character. 
For V = 2.2 the coherent QP band becomes flatter and an 
incoherent band appears in the energy range u ~ 0.5t — 2.5t. 
The QP peak has been cut for clarity purposes. 



contains information of both, the density of states of the 
interacting boson and its coupling with fcrmions but, in a 
mixed form. For instance, the only way that a 2 F follows 
exactly the same shape of the boson density of states oc- 
curs if the coupling does not depend on u> or momentum. 
This condition is not satisfied in our case which can be 
seen by looking at the expression (fTT)]) for Ira's. On the 
other hand, only in the first term of eq. (fT0|) the charge- 
charge correlation is explicitly present. The other terms 
which contain D\r and D\\ are proper of our method, 
and they are due to the nondouble occupancy constraint. 

With increasing V, the features marked with solid ar- 
rows in Fig.[5]become soft and, at the same time, spectral 
weight appears at lower energies (dashed arrows). This 
behavior is clearly correlated with that depicted in Fig. 
[TUl which shows the soft charge dynamics discussed in 
Sec. EH (see Fig.©. 

It is well known2. 4 that 

-^=L^M. (i 5) 

Therefore, the presence of low energy excitations in a 2 F 
leads to an increase of dRe T,/du> and hence, to a decrease 
of the QP weight Z. Notice that near V c , a 2 F(uj) shows, 
besides the low energy features, structure at high energy 
in agreement with the behavior presented by charge cor- 
relations (Fig. Hand Fig. [lOj). The fact that not all the 
spectral weight is concentrated at low energy is the cause 
for the slow decrease of Z. The inset of Fig. [5] shows 
that Z decreases somewhat faster for x = 2/3 than for 
x = 1/3 which is consistent with the fact that softening 
of the charge collective mode is more clear for x = 2/3 
(see Figs. and 0}. 



To conclude, we think that we have given clear argu- 
ments which show that the charge collective soft modes 
near the CDW are responsible for the incoherent motion 
of the carriers near charge ordering. In addition, the 
charge soft modes reach lu = at V c which means that 
—dRe Yj/doj — > oo (eq. [T5|) and then Z — > as indicated 
in Fig. M 



E. E-k structures 

When discussing the spectral functions we have shown 
the presence of coherent QP peaks and incoherent struc- 
tures. In Fig.[TT]we show the energy position of the main 
visible structures for V — and V^ = 2.2ata; = l/3 along 
the main directions of the Brillouin zone (inset in Fig. 
[TTj) . Dashed lines mark the spectral functions at kp. For 
V = (panel (a)), we can see that the QP band (near 
the FS) and a dispersing incoherent band at lu ~ — 3t 
dominate the spectra. 

When V increases to V = 2.2 (panel (b)) the QP co- 
herent band becomes less dispersing consistently with the 
reduction of the QP weight. The incoherent band at neg- 
ative u> becomes also less dispersing. Interestingly, an 
incoherent band is formed in the region ~ 0.5i — 2.5<. 
Notice that no kink is showed by the coherent QP band. 
This is because the kink, if it exist, should be located at 
energies of the order of the interacting boson (which in 
this case is related to the soft mode) and, at V = 2.2, it 
is still located at energies of ~ 0.2t (Fig. 2]) which is of 
the order of the QP bottom energy. 
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V. COMPARISON WITH COBALTATES 

In previous sections we have discused one-particle spec- 
tral predictions of the large-TV approach for the t-J-V 
model on the triangular lattice without making emphasis 
on any particular system. In this section we will discuss 
previous results and their possible contact with cobal- 
tates. First principle calculations in cobaltates^ predict, 
besides a large hole Fermi surface around T, the existence 
of six small pockets near the corners of the Brillouin zone. 
Until now, all ARPES experiments both, in hydrate d 27 ' 28 
and unhydrated cobaltates^^^^i^ 2 -^^ do not show 
the presence of these small pockets. The absence of pock- 
ets in ARPES can be understood by the renormalization 
of the bands due to strong electronic correlations 35 . In 
fact, a reduction of the bandwidth and large electronic ef- 
fective masses are observed in ARPES 31 ' 32 . In addition, 
in a very recent work 3 ^ it was pointed out that cobaltates 
are close to a Mott insulator in the limit x — > 

Following these ideas, a one-band t-J model on the tri- 
angular lattice, as proposed in the present paper, seems 
to be appropriated for studying cobaltates. For this pur- 
pose, c\ a and Ci<j in eq. dXJ) must be associated with the 
creation and destruction of holes respectively, as in Refs. 

Then, in this context, the doping x is the elec- 
tron doping away from half-filling. Considering the bare 
hoping t ~ l50meV— and Udd ~ 3.0 — 5.5eU— we obtain 
J/t ~ 0.1—0.2, similarly to the value used in previous sec- 
tions. In addition, we will mainly focus the comparison 
with our results for x — 1/3, which is close to the doping 
where the maximum superconducting critical tempera- 
ture T c is found in cobaltates. For x < 0.5 our results 
predict a Pauli paramagnet as observed in cobaltates. On 
the other hand, as cobaltates exhibit a Curie- Weiss be- 
havior in the region x > 0.5, we will not relate our results 
in a; = 2/3 with experiments. However, it is important to 
mention that it is not clear if this Curie- Weiss behavior 
is due to the fact that electronic correlations are more 
important for x > 0.5 or other effects, as for instance Na 
order, are predominant (see Ref. [401] and references their 
in). 

According to results in previous sections we have two 
factors which reduce the Fermi velocities from the bare 
one. The first one occurs already at mean field level 
where the effective hopping is (x/2)t instead of t. The 
second one, which is obtained only after evaluating fluc- 
tuations beyond mean field is the quasiparticle weight 
Z. Therefore, the Fermi velocity is vf ~ ( x /2) Z Vp DA 
where Vp DA is the bare Fermi velocity extracted from 
LDA calculations. Using the parameters for cobaltates 
and our results for x = 1/3 up to 0.9U C we obtain 
vf ~ 0.21eVA This value is close to the value vp — 
(0.30±0.05)eUl found in Refs. This small value 

was considered 32 as an indication for the presence of 
strong correlations and, the ratio between T c and the 
Fermi velocity was found to be close to the same ratio 
for cuprates. In agreement with these experiments, our 
calculation shows also a rather isotropic Fermi velocity 



over the BZ (see Fig. fTTj) . Another interesting and quali- 
tative agreement between experiments and our finding is 
the following: in Ref. [HJ it was found that the scatter- 
ing rate behaves as ~ lo 2 near the Fermi surface which 
agrees with results in panel (a) of Fig. [7l 

In our previous papers^, for describing superconduc- 
tivity in cobaltates we proposed, like other work s 3 ' 4 ' 37 , 
that the system is close to V3 x \/3 charge or- 
der. The charge order scenario is supported by some 
experiments 7 -^. In addition, ARPES experiments 3 -^ 
show a Fermi surface topology, very close to our FS, 
which favors charge instabilities. We have proposed 
that the interplay between electronic correlations and 
phonons are relevant for describing superconductiv- 
ity with triplet NNN-/ pairing as predicted by some 
experiment s 41 ' 42 ' 43 ' 44 ' 45 . For this scenario we need that, 
with increasing water content, the system comes closer to 
charge ordering. Recent ARPES experiment s 27 ' 28 have 
shown that the Fermi velocity do not vary very much 
with the water content. These experiments seem to be in 
contradiction with our result of Fig. [TTJ where it is pos- 
sible to see that very close to the charge instability the 
band dispersion becomes flat. However, we think that 
it is possible that for the water content reached in the 
experiments the system is still far from the static charge 
order state i.e., in the regime where Z does not change 
very much. This could explain the reason why NMR 46 
does not show the charge order because it is fluctuating 
faster than the time-scale accessible by NMR. We remem- 
ber that for obtaining a robust description of supercon- 
ductivity we do not need to place the system very close 
to the charge order being enough a V ~ 80 — 90%V^&. 
However, for V > 0.9V C , Z is strongly reduced and two 
competing effects can be expected. From one side, the 
decreasing of Z may cause an increasing of the electronic 
density of states and then, a reinforcement of supercon- 
ductivity. On the other side, when Z decreases the quasi- 
particle losses coherence thus, superconductivity may be 
diminished. A careful study of this competition requires 
the calculation of superconducting pairing in 0(1/N 2 ). 

It is interesting to compare our results with a very re- 
cent ARPES experiment on a series of Bi cobaltates^ 7 -. 
These compounds contain the same triangular Co planes 
as Na cobaltates. Besides a flat QP band with similar 
Fermi velocity as in Na cobaltates, incoherent features 
were observed at high energy in the Bi cobaltates. In 
that paper it was discused that these incoherent features 
are created with weight transfered from the QP band, a 
mechanism similar to that discussed in Sec. IIVI (notice 
the analogy, for instance, between Fig. 2 of Ref. [47] and 
our Fig. ITT)) . On the other hand, the center of the inco- 
herent band in Fig. 2 of Ref. [47| is at ~ —0.5 eV, close 
to the center of our incoherent band (Fig. [TTj) which, us- 
ing t — 150 meV is about —0.45 eV. Similar high energy 
features were recently observed in cuprates (see Ref. [47| 
and references their in) and, the present approach was 
used 48 for discussing those features occuring in the over- 
doped side where, as mentioned before, our method is 
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expected to be reliable. 

Optical conductivity experiment a 49 i 50 show a different 
behavior to usual metals where, besides a Drude peak at 
low energies, a broad absorption centered around ~ 250 
meV is observed and interpreted as a pseudogap 5 ^ be- 
havior. Similar features are observed in organic materi- 
als and discussed in term of the charge order proximity 51 . 
From our results in Fig.[5]we expect that, near the charge 
order, transitions from the Drude peak to the broad 
structure formed between ~ 75 — 375 meV (where we 
used t = 150 meV) are possible leading to the optical 
absorption as in the experiments. 

Recently, optical conductivity experiments in organic 
materials^ have also shown an anomalous increa sing of 
the effective mass with temperature. In Ref. [22|, a 
reentrant behavior in the T-V phase diagram, as the 
one discussed in our Fig. [5j was found to be responsi- 
ble for the increase in effective mass. If for cobaltates 
0.7V C < V < V c , we think that an increase in effective 
mass with temperature can also be expected in these ma- 
terials. From Fig. [51 depending on the value of V, the 
effective mass increase may occurs for T < 800 K. It will 
be interesting to perform this kind of experiments. 

VI. CONCLUSION AND DISCUSSIONS 

In this paper, the t-J-V model at finite density was 
studied on the triangular lattice. V represents the 
Coulomb interaction between nearest neighbors and it 
is the parameter responsible for triggering a charge or- 
der phase for V > V c . V c is the critical Coulomb re- 
pulsion whose numerical value depends on doping (Fig. 
[5]). Studying charge correlation functions, it was shown 
that near V c charge dynamics becomes soft and, at 
V = V C , charge modes collapse to zero frequency freezing 
a a/3 x y/3 CDW phase (Figs. [3] and 0J). 

By evaluating fluctuations beyond the mean field level, 
one-particle spectral properties were computed. Far away 
from charge ordering, i.e. for zero or small V, the QP 



weight Z is reduced from one due to pure t-J model cor- 
relations effects. With increasing V , Z decreases slowly 
until V ~ 90% of V c - Beyond this value Z decreases 
faster approaching zero at V c (Fig. [5]). Therefore, near 
V c the electronic dynamics becomes very incoherent. For 
V < V c , the scattering rate behaves as ImH ~ ui 2 (Fig. 
Eta)) which is characteristic for a Fermi liquid behavior. 
In addition, when the system approaches charge ordering, 
Im X accumulates spectral weight at low energy. Due to 
this fact, the slope of Re S at u> — increases (Fig. [7](b)) 
leading to the decrease of Z discussed above. 

From Im X and Re X the spectral function A(k, to) was 
calculated. With increasing V, approaching V c , the most 
important changes are present at energies close but larger 
than the Fermi energy (w = 0). While the weight of the 
QP peak decreases (Z decreases), spectral weight appears 
at small positive energy in the range ~ 0.5i — 2.5i (Fig. 
IS). 

The calculation of a 2 F(ui) allows to identify the exci- 
tations responsible for the self-energy effects. Comput- 
ing a 2 F(u>) as a function of V approaching V c , it was 
shown that its behavior can be, without ambiguity, cor- 
related with charge spectra (Figs. [91 and [TO)) . Therefore, 
the charge soft modes, responsible for the charge instabil- 
ity, lead to an increase of a 2 F at small frequencies with 
the corresponding decrease of Z. 

In Sec.fVl we mainly confronted our results with recent 
ARPES experiments in cobaltates. We found our results 
to be in agreement with experiments presented in Refs. 

Finally, we would like to point out that while our re- 
sults were obtained for a purely 2D system, we expect 
only quantitative differences in our results when passing 
to a corresponding anisotropic 3D system, since we are 
dealing here with a spontaneous symmetry breaking of 
a discrete symmetry, as corresponds to a commensurate 
CDW. Certainly critical exponents may change, but they 
are beyond the scope of our treatment. 
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